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I.  INTRODUCTION 


The  NASTRAN1,2  (NASA  Structural  Analysis)  finite  element  program  is 
a  general-purpose,  3TgitaT7  structural-analysis  computer  program  which 
uses  the  finite  element  method  to  analyze  the  behavior  of  elastic 
structures.  This  computer  program  is  in  use  at  the  Ballistic  Research 
Laboratory  (BRL)  on  the  UNIVAC  1108.  Among  the  problems  requiring 
solution  at  the  BRL  are  many  that  Involve  rectangular  plates  and 
membranes.  Some  of  the  problems  require  the  generation  of  a  two- 
dimensional  grid  while  others  may  require  a  three-dimensional  grid. 
Some  problems  require  using  the  CQUAD2  plate  elements  to  obtain  a 
solution.  If  the  plate  is  large,  a  great  number  of  grid  points  and 
CQUAD2  elements  may  be  needed  to  adequately  define  the  structure.  The 
result  is  a  considerable  amount  of  figuring,  writing  and  then  punching 
of  cards,  all  of  which  is  time-consuming  and  error-prone  for  large 
systems.  Such  a  large  problem  resulted  in  two  programs,  GRID  and  CQUAD, 
being  designed  to  generate  the  required  two-dimensional  grid  and  the 
CQUAD2  connecting  elements,  respectively,  and  to  punch  the  cards 
required  for  the  grid  and  element  input.  The  GRID  program  was 
subsequently  modified  to  be  able  to  generate  a  three-dimensional  grid. 

II.  DISCUSSION 

The  GRID  program  is  listed  in  Figure  1.  The  input  to  the  program 
consists  of  two  cards.  The  first  card  gives  A,  B,  and  C,  the  x,  y,  and 
z  spacings  of  the  grid,  respectively,  in  FORMAT  (3F10.4).  The  second 
card  gives  X9,Y9,  and  29,  the  x,  y,  and  z  dimensions  of  the  grid,  in  the 
same  format.  In  the  program,  the  two  formats  are  given  different 
numbers  to  leave  open  the  option  of  changing  the  format  of  one  of  the 
sets  in  a  simple  manner.  As  may  be  noted  the  fields  3,  7,  and  8  are 
left  blank  since  they  can  be  specified  on  a  GRDSET  card,  or  permanent 
single-point  constraints  associated  with  specific  grid  points  could  be 
punched  in  later.  Fields  3  and  7  refer,  respectively,  to  the 
identification  number  of  the  coordinate  system  in  which  the  grid  point 
is  defined  and  to  the  identification  number  of  the  coordinate  system  in 
which  the  displacements  are  measured  at  the  grid  point;  field  8 
contains  the  single-point  constraints  that  are  to  restrict  the  motion  of 
the  grid  point.  If  a  two-dimensional  grid  is  desired  then  29  is  left 
blank,  and  the  program  will  then  generate  a  two-dimensional  grid  and 
punch  out  the  appropriate  GRID  cards.  If  the  uniform  spacing  does  not 
cover  the  range  of  the  grid  exactly,  then  the  remainder  is  placed  at  the 
end;  i.e.,  2.75  units  if  a  1-unit  spacing  is  specified  would  result  in 
two  1-unit  spacings  and  the  last  one  with  a  0.75  spacing.  This  is  true 
for  each  of  the  dimensions,  x,  y,  and  z.  Figure  2  gives  the  resultant 
grid  listing  where  a  three-dimensional  grid  had  X9,  Y9,  29,  A,  B,  and  C 


^he  NASTRAN  Theoretical  Manual,  NASA  SP-222(01),  December  1972 
2 The  NASTRAN  User's  Manual,  NASA  SP-22(01),  May  1973 


7 


ooooooooooo 


PROGRAM  GRID (INPUT*0UTPUT*PUNCH*TAPE5*INPUT» 
1TAPE6*0UTPUT*TAPE7*PUNCH) 

X9,Y9fZ9  ARE  THE  X ,Y ,  AND  Z  DIMENSIONS  OF 
THE  PLATE.  IF  A  TWO  DIMENSIONAL  GRID  IS 
DESIRED  *  PUT  A  ZERO  IN  FOR  Z9#  A»B»C  ARE 
THE  UNIFORM  DIMENSIONS  OF  THE  ELEMENTS  IN 
THE  X.YfZ  DIRECTIONS*  RESPECTIVELY.  IF  THE 
GRID  DOES  NOT  COME  OUT  EVEN  IN  ANY  DIRECTION 
THEN  THE  REMAINDER  IS  PUT  IN  THE  LAST 
INCREMENT  OF  THAT  DIRECTION.  IT  IS  THERE¬ 
FORE  PREFERABLE  TO  SELECT  A*B«  AND  C  TO 
MAKE  THE  GRID  SPACINGS  UNIFORM  OVER  THE 
ENTIRE  GRID. 

INTEGER  D*E*F*S 
READ (5*50) A«B«C 

50  FORMAT (3F10. 4) 

WRITE  (6  *5  0  )  A  *B  *C 
X«=  0 

Y  =  0 
Z*0 

READ(5»51) X9»Y9*Z9 

51  FORMAT (3F10. 4) 

WR ITE ( 6*  51  )X9* Y9»Z9 
S=1 

120  FORMAT (8H  GRID  *  1 8* 8X « 3F8.3) 

L=INT (X9/A+1 ) 

M*  INT ( Y9/B  *1 ) 

N=1 

D*M 

IF (Z9. EO.O) GO  TO  168 
NSIN  T ( Z9/C  + 1 ) 

168  CONTINUE 

DO  200  K*1»N 

X=0 

Y*0 

169  CONTINUE 

170  DO  400  Jsl,M 

179  X«0 

180  CONTINUE 

DO  500  I=1*L 
WRITE  (6* 120 ) S»X» Y  »Z 
WRITE  (7  *  120) S*X*Y«Z 
S*S*1 
XSX  + A 


Figure  1.  Listing  of  GRID  Program. 
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500  CONTINUE 
X*X-A 

IF  (X.NE.X9)G0  TO  503 
GO  TO  504 

503  X*X9 

WRITE  <6tl20)StX,Y*Z 
WRITE  (7t 120) S*X* Y»Z 
S*S*1 

504  Y*Y+B 
400  CONTINUE 

Y  =  Y-B 

IF  (Y.EQ.Y9) GO  TO  190 
CONTINUE 
402  Y*Y9 

GO^TO  170 
190  Z=Z*C 
M=D 

200  CONTINUE 
Z=Z-C 

IF(Z.EQ.Z9)G0  TO  700 

Z*Z9 

M=D 

N=1 

GO  TO  168 

700  STOP 

701  END 


Figure  1  Continued.  Listing  of  GRID  Program. 
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GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GPID 

GRID 

GRID 

GRID 

GRID 

GPID 

GRID 

GRID 

GPID 

GRID 

GRID 

GRID 

GRID 

GPID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GPID 

GPID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 

GRID 


1 

0.000 

0.000 

0.000 

2 

1.000 

0.000 

0.000 

3 

2.000 

0.000 

0.000 

4 

2.750 

0,000 

0.000 

5 

0.000 

1.000 

0.000 

6 

1.000 

1.000 

0.000 

7 

2.000 

1.000 

0.000 

8 

2.750 

1.000 

0.000 

9 

0.000 

2.000 

0.000 

10 

1.000 

2.000 

0.000 

11 

2.000 

2.000 

0.000 

12 

2.750 

2.000 

0.000 

13 

0.000 

2.500 

0.000 

14 

1.000 

2.500 

0.000 

15 

2.000 

2.500 

0.000 

16 

2.750 

2.500 

0.000 

17 

0.000 

0.000 

1.000 

18 

1.000 

0.000 

1.000 

19 

2.000 

0.000 

1.000 

20 

2.750 

0.000 

1.000 

21 

0.000 

1.000 

1.000 

22 

1.000 

1.000 

1.000 

23 

2.000 

1.000 

1.000 

24 

2.750 

1.000 

1.000 

25 

0.000 

2.000 

1.000 

26 

1.000 

2.000 

1.000 

27 

2.000 

2.000 

l.ooo 

28 

2.750 

2.000 

1.000 

29 

0.000 

2.500 

1.000 

30 

1.000 

2.500 

1.000 

31 

2.000 

2.500 

1.000 

32 

2.750 

2.500 

1.000 

33 

0.000 

0.000 

1.750 

34 

1.000 

0.000 

1.750 

35 

2.000 

0.000 

1.750 

36 

2.750 

0.000 

1.750 

37 

0.000 

1.000 

1.750 

38 

1.000 

1.000 

1.750 

39 

2.000 

1.000 

1.750 

40 

2.750 

1.000 

1.750 

41 

0.000 

2.000 

1.750 

42 

1.000 

2.000 

1.750 

43 

2.000 

2.000 

1.750 

44 

2.750 

2.000 

1.750 

45 

0.000 

2.500 

1.750 

46 

1.000 

2.500 

1.750 

47 

2.000 

2.500 

1.750 

48 

2.750 

2.500 

1.750 

Figure  2.  Three-Dimensional  Grid  Generated  by  GRID  Program. 
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specified  as  2.75,  2.5,  1.75,  1.0,  1.0,  and  1.0.  Figure  3  gives  the 
resultant  grid  listing  where  a  two-dimensional  grid  had  X9,  Y9,  2)*  A, 
B,  and  C  specified  as  3.75,  7.75,  blank,  1.0,  1.0,  and  blank. 

Figure  4  gives  the  listing  of  the  CQUAD  program.  The  Input  to  the 
program  consists  of  3  Integers  In  FORMAT  (3110),  giving  the  total  number 
of  grid  points  In  the  mesh,  the  number  of  grid  points  In  the  x  direction 
and  the  Identification  number  of  the  PQUAD2  property  card.  Field  8,  the 
material  property  orientation  angle, Is  left  blank.  The  CQUAD  program 
creates  the  grid  connections  In  a  counter-clockwise  direction  around  the 
perimeter  of  each  element  and  numbers  the  elements  from  one  to  the  total 
number  of  CQUAD2  elements.  Figure  5  gives  a  listing  of  the  CQUAD2  cards 
generated  using  the  grid  from  Figure  3. 

III.  CONCLUSION 

The  two  programs  have  been  tested  on  the  CYBER  7600  for  two-  and 
three-dimensional  grids  and  for  the  dimensions  being  exact  numbers  of 
Intervals  and  for  grids  having  dimensions  not  an  exact  number  of 
Intervals.  The  CQUAD  program  was  also  tested  on  the  CYBER  7600.  Use  of 
these  programs  will  reduce  both  time,  labor,  and  the  probability  of 
error. 
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3 
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GRID 
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GRID 

5 
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GRID 
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GRID 

8 
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9 

3.000 

1.000 

0.000 
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11 

0.000 

2.000 
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GRID 

12 
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GRID 
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GRID 

14 

3.000 

2.000 

0.000 

GRID 

15 

3.750 

2.000 

0.000 

GRID 

16 

0.000 

3.000 

0.000 

GRID 

17 

1.000 

3.000 

0.000 

GRID 

18 
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3.000 

0.000 

GRID 

19 

3.000 

3.000 

0.000 

GRID 

20 

3.750 

3.000 

0.000 

GRID 

21 

0.000 

4.000 

0.000 

GRID 

22 

1.000 

4.000 

0.000 

GRID 

23 

2.000 

4.000 

0.000 

GRID 

24 

3.000 

4.000 

0.000 

GRID 

25 

3.750 

4.000 

0.000 

GRID 

26 

0.000 

5.000 

0.000 
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27 
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5.000 

0.000 

GRID 

28 
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5.000 

0.000 

GRID 

29 

3.000 

5.000 

0.000 

GRID 

30 

3.750 

5.000 

0.000 

GRID 

31 

0.000 

6.000 

0.000 

GRID 

32 

1.000 

6.000 

0.000 

GRID 

33 

2.000 

6.000 

0.000 

GRID 

34 

3.000 

6.000 

0.000 

GRID 

35 

3.750 

6.000 

0.000 

GRID 

36 

0.000 

7.000 

0.000 

GRID 

37 

1.000 

7.000 

0.000 

GRID 

38 

2.000 

7.000 

0.000 

GRID 

39 

3.000 

7.000 

0.000 

GRID 

40 

3.750 

7.000 

0.000 

GRID 

41 

0.000 

7.750 

0.000 

GRID 

42 

1.000 
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0.000 

GRID 

43 
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0.000 
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44 
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0.000 
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45 

3.750 

7.750 

0.000 

Figure  3.  Two-Dimensional  Grid  Generated  by  GRID  Program. 
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o  o  o  o 


G9  IS  THE  TOTAL  NUMBER  OF  GRID  POINTS  AND 
G8  IS  THE  NUMBER  OF  GRID  POINTS  IN  THE  X 
DIRECTION.  P  IS  THE  PROPERTY  ID  FOR  THE 
CQUAD2  ELEMENT. 

INTEGER  S.P.0.Q.G9.G8 
S  =  0 
M  =  0 

READ(5»5Q)G9»G8»P 
WRITE (6.50JG9.G8.P 
70  INDEX=G8-1 

DO  200  1=1. INDEX 

S  =  S*1 

M=M*1 

N=M*  1 

0=N+G8 

0  =  0-1 

WRITE (6.20) S.P.M.N.O.Q 
WRITE (7.21) S.P.M.N.O.Q 
200  CONTINUE 

IF  (0 . EQ .G9 ) GO  TO  180 
M=M+1 
GO  TO  70 
180  STOP 

21  FORMAT (8HCQUAD2  .618) 

20  FORMAT ( 9H  CQUAD2  .618) 

50  FORMAT (3110) 

END 


'Figure  4.  Listing  of  CQUAD  Program. 
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Figure  5.  CQUAD2  Elements  Generated  by  CQUAD  Program  for  Grid  of 
Figure  3. 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number _ _ ___ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _  _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


